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NAVAL ARCHITECTURE 


Investigations of influence 
of screw propeller operation on water flow 
around stern part of ship hull 


Tadeusz Koronowicz 
Zbigniew Krzemianowski 
Institute of Fluid-Flow Machinery, 
Polish Academy of Sciences 

in Gdansk 


ABSTRACT 


This paper presents results of measurements of velocity field in before- the - propeller flow in presence of 

a ship model hull of two configurations, as well as comparative calculations of velocity field on a full-scale 

ship. Analysis of the research results showed that input data to Biot-Savart formula should be modified in 
the case of calculations of propeller-induced velocities on ship hull surface. 


Keywords : ship hydromechanics, propeller-induced velocities, Biot-Savart equation. 


INTRODUCTION 


In the 1990s in Ship Propeller Division, Institute of Fluid- 
-Flow Machinery, Polish Academy of Sciences (IMP PAN), the 
computer model basin PANSHIP was elaborated. It has been 
aimed at simulation of ship hull model tests in ship model basin 
as well as relevant calculations for full-scale ship hull. 

The computer model basin is a computer software system 
consisted of a dozen or so mutually cooperating programs 
[1+6]. Crucial elements of the system are the programs capable 
of taking into account the influence of screw propeller operation 
on flow around ship hull. The software contains the programs 
with the use of which a change of hull resistance resulting 
from propeller suction action can be determined, and those by 
which the influence of propeller operation on velocity field in 
behind-the-hull flow can be taken into account. 

The initial calculations performed by using the PANSHIP 
software have yielded generally correct results with the excep- 
tion of one element : changes of hull resistance resulting from 
propeller operation. In ship theory such change is expressed 

in the form of the so-called thrust deduction t : 


t=(R,-R,V/T=(T-RYT 


where : 
R, — resistance of hull without propeller 
T — propeller thrust 
R, — resistance of hull with operating propeller (identified 
with propeller thrust) 


The quantity t is usually determined during every ship model 
propulsion tests in model basin. 


As such tests have been performed every year for many 
ship models, a very rich collection of experimental data in this 


domain has been gathered. Basing on them one can unambigu- 
ously state : the more full form of a ship the greater value of 
its thrust deduction t. 

In the preliminary version of the computer model basin 
in question, for hulls of more full forms, greater and greater 
differences between calculated values of thrust deduction and 
those experimentally determined for the same hulls, were 
obtained (Fig.1). Due to prior research on flow around propel- 
ler it was possible to diagnose that the velocities induced by 
whirls representing the propeller, determined by means of 
the Biot-Savart formula, obtained erroneous values on the 
hull surface. In Fig.2 it can be observed that the more full form 
of a hull the smaller values of the induced velocities calculated 
from the original Biot-Savart formula, therefore the calculated 
pressures (under-pressures) on the hull surface take also smaller 
values. 


0.3 t 
0.25 

0.2 | 

Test Calculation with a modified 
0.15 Biot-Savart equation =| 
0.1 
Original calculation 

0.05 

0 ô 

0.4 0.45 05 0.55 06 065 07 0.75 08 0.85 0.9 


Fig. 1. A simplified diagram of the relation between the thrust deduction 
and hull block coefficient (in reality the relation is more complex 
as it depends first of all on fullness of stern part of hull) . 


By analyzing the diagrams presented in Fig.2 and 3 
it can be explained why such results have been obtained. 
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0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 


Fig. 2. Distribution of axial component of velocity in before-the-propeller 
flow, obtained from the original Biot-Savart formula, shown on the backgro- 
und of frame sections of hulls having various block coefficients . 


If the velocities induced by the whirl system which repre- 
sents propeller blades and propeller race, are calculated from 
the original Biot-Savart formula then the velocity distribution 
starts at the hull plane of symmetry (the point C’ in Fig.3). 
The similar velocity distribution is presented in Fig.2, where 
simultaneously the frame sections of 3 ship hulls of different 
values of the block coefficient 6 are shown. It can be observed 
that the greater fullness of the hull the smaller obtained values 
of velocities induced on its surface. 


Fig. 3. Schematic presentation of the modification 
of data input to the Biot-Savart formula. 


In the up-to-date version of the PANSHIP, was implemented 
a new method of calculation of induced velocities by means 
of Biot-Savart formula, (called the engineering method). The 
propeller-induced velocities were calculated in the point C’ 
(Fig.3), but they were considered as the velocities calculated 
on the hull (the point C on the hull). It means that zero-value of 
the coordinate perpendicular to the hull plane of symmetry was 
put in the Biot-Savart formula. The calculation results appeared 
significantly better. Values of the thrust deduction obtained from 
calculations and those from experiments became more and 
more similar to each other. Obviously the described method of 
determination of induced velocities is only approximate, how- 
ever, as revealed from practice, it yielded satisfactory results 
in engineering applications without any special modifications 
of the software. 

For many years the so-modified computer software PAN- 
SHIP has been in use, and the hypothesis associated with 
the modification of input data for Biot-Savart formula was 
confirmed by comparing calculation results with experimental 
ones. However it was necessary to test the hypothesis by means 
of direct measurements of the field of the propeller - induced 
velocities around the hull. Such a verification is the subject of 
the presented work. 
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The model tests were performed at the Ship Hydromecha- 
nics Centre of CTO [10]. They consisted in measuring the velo- 
city field around stern part of ship both without any propeller 
and with operating propeller. 

The investigations were conducted on the ship model 

having its main particulars as follows : 


Length b.p. —-— 6515m 
Breadth — 0.977m 
Draught — 0.376m 
Model scale -— A=33 


The hull frame sections are shown in Fig.4. 


Bow 


Stern 
Fig. 4. Image of the panels projected on the model frame sections . 


The applied measuring instrument (the measuring sounder 
fitted with the single five-hole spherical head PKN(5+4)/8/1) 
made it possible to measure the velocity components Vx, Vy, 
Vz in the hull-fixed rectangular coordinate frame : 


> x-— component: along ship axis of sym- 
metry and hull motion direction 

> y- component : perpendicular to the 
hull plane of symmetry 

> z— component : perpendicular to the 
hull water plane. 


The measurement space was located 
at the port side of the hull. 


Distance of propeller working plane from 
aft perpendicular — Xp = 124 mm 


Distance of propeller axis from 
plane of symmetry — Yp = 0.0 mm 


Distance of propeller axis from 
base plane — Zp = 109.1 mm. 


The measurement plane was located AX = 157 mm 
fore from the propeller working plane. 


EXPERIMENTAL TESTS 


The measurements were performed at one value of the ship 
model velocity V,, = 1.75 m/s and four values of rotational 
speed of the propeller model. The first value of rotational 
speed was determined for zero-value of propeller thrust. It 
was assumed that this was the rotational speed at which pro- 
peller-induced velocities on the hull were of negligibly small 
values, hence the velocity measurements could be considered 
equivalent to the tests on the hull without propeller. The value 
of n = 7.3 rps resulted from the tests (during the tests values 
of both propeller - induced thrust and torque as well as hull 
resistance were measured). 

The next three values of rotational speed were so selected 
as to obtain only significantly large values of propeller - model- 
-induced velocities. With taking into consideration the working 


range of the measuring dynamometer the following three values 
of rotational speed were selected : 


n,=25 l/s ; n,= 30 1/s 
n,=35 Tis, 


At the obtained values of rotational speeds the values of 
propeller thrust were many times greater than that of hull re- 
sistance at the speed V,,= 1.75 m/s. 


The measurements were conducted along two measurement 
lines perpendicular to the longitudinal plane of symmetry of 
the hull, XZ, placed by AX = 157 mm apart from the propel- 
ler working plane. The first line is placed at the height of 
propeller rotation axis, the other - 50 mm above the men- 
tioned axis. 


The first measurement point was selected as close to the 
hull surface as possible, and the successive points were placed 
at every 20 mm up to the distance assumed negligible from the 
point of view of propeller — induced velocities. 


The selected measurement results are presented in Fig. 5+8 
whereas the complete set of them - in CTO’s report [10], and 
their graphical representation - in the IMP PAN report [11]. 
Values of the velocity components Vx, Vy, Vz and of the total 
velocity Vc can be found there. During all the tests the ship 
model speed V,, was kept equal to 1.75 m/s. 

In Fig.5 are presented the measurement results at the 
rotational speed n = 7.3 1/s corresponding with zero-value 
of propeller model thrust. Hence it can be assumed that the 
velocities shown in Fig. 5 correspond with those around the 
hull without propeller. They have been taken as the reference 
point for determining the velocities induced by working pro- 
peller model. 

2 
1.8 
1.6 
1.4 
1.2 

1 
0.8 
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Fig. 5. Velocity components along the measurement line located at the 
height of the propeller axis, for n = 7.3 1/s (thrust of zero-value) . 


V [m/s] 


In Fig.6, 7 and 8 are presented results of the measurements at 
higher rotational speeds of propeller model, for which induced 
velocities should already show significant values. In Fig.9 it 
can be observed in which way values of the axial component 
(marked x) change along with rotational speed changing. 
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Fig. 6. Velocity components along the measurement line located 
at the height of the propeller axis, for n = 25 I/s . 
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Fig. 7. Velocity components along the measurement line located 
at the height of the propeller axis, for n = 30 I/s . 
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Vm = 1.75 m/s — 
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Fig. 8. Velocity components along the measurement line located 
at the height of the propeller axis, for n = 35 I/s . 


V, [m/s] 


n= 35.0 1/s 
1.47—— n = 30.0 I/s 


1.2 ->——- n= 25.0 1/s 


Vm = 1.75 m/s — 


0.8 +n =_ 7.3 I/s 


0 5 10 15 20 ylem] 25 


Fig. 9. Modules of velocities along the measurement line located 
at the height of the propeller axis, for various values 
of propeller rotational speed . 


In Fig.10 are presented the differences between the veloci- 
ties obtained at high values of rotational speed and the velocity 
corresponding with the thrust of zero-value. They should cor- 
respond with the propeller-induced velocities but the character 
of the changes indicates that the influence of viscosity on the 
velocity distribution is significant (induced velocities make 
velocity distribution in the boundary layer changing). 


0.7 


4 T 
Vi [m/s] y= 35.0 1/s Vic 
0.6 + + = 
n= 30.0 1/s 

0.5 + 
0.4 n = 25.0 1/5 œ Vm=1.75 m/s__| 
0.3 
0.2 
0.1 

0 

0 5 10 15 20 y[em] 25 


Fig. 10. Axial component and module of propeller-induced velocity along 
the measurement line located at the height of the propeller shaft axis . 
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Similar measurements were performed at the measurement 
line located in the same plane but at the height above the pro- 
peller rotation axis by 50 mm. For the measurements only the 
final diagram of the propeller induced velocities is presented. 


0.7 


. i f 

Vi [m/s] Pa n=35.0 I/s 

0.6 Ñ 
n = 30.0 1/s 
0.5 = 
Vix 
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0.3 
0.2 
0.1 
0 
0 5 10 15 20 25 30 35 


Fig. 11. Axial component and module of propeller-induced velocity 
along the measurement line located at the height 
above the propeller shaft axis by 50 mm . 


Analyzing the above presented results of the investigations, 
especially those of Fig. 10 and 11, one can state that the share 
of propeller-induced velocities in the total velocity is significant 
(intentionally the values of propeller rotational speed conside- 
rably exceeded the own propulsion point of the model, which, 
for the velocity V=1.75 m/s, approximately corresponded with 
the rotational speed of 13 rps). 

To confirm hull influence on calculation results of propeller- 
-induced velocities it should be necessary to perform measure- 
ments in two frame cross-sections located nearby to each other 
but having significantly different transverse offsets (breadth). 
Unfortunately, for many years the ship models designed and 
tested have been characterized by slender stern forms. For 
this reason the second cross-section was chosen beyond the 
stern. In order to maintain the distance between the working 
plane and measurement cross-section the way of fastening the 
propeller shaft was changed. The stern part of the hull was 
modified by extending the stern tube in such a way as to get 
the measurement cross-section placed beyond the stern and the 
propeller model placed at the same distance as in the case of 
basic tests (Fig.12). The tests on such hull version were called 
the tests on the Ship 2. 


a 


Z 100.50) 


V7 
Fig. 12. Location of measurement lines on the ship model 
with the modified stern. PP - Base plane, PR - Aft perpendicular . 


The measurements were performed at the same ship model 
speed V,, = 1.75 m/s and four values of rotational speed of 
propeller model, in the same way as in the first cycle of inves- 
tigations. Their results for the sounder located at the height of 
the shaft axis are presented in Fig.13. (The comprehensive set 
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of the results from the measurements and calculations can be 
found in the CTO report [10] and IMP PAN report [11]). 
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Fig. 13. Axial component and module of propeller-induced velocity 
measured at the height of the propeller shaft axis of the ship model 
with the modified stern, at the distance of 157 mm before the propeller . 


Comparing, with each other, the measurement results for 
both versions of propeller fastening (Fig. 14) one can state that 
the curves are mutually shifted. The difference is approximately 
equal to the difference of hull breadth and propeller shaft in the 
places where the measurements have been performed in both 
versions of the tests. 


0.7 T 
Vic [m/s] n=35.0 1/s 
0.6 Vm = 1.75 m/s 
n= 30.0 1/s tests on the level 
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n= 30.5 I/s 
n=25.5 1/s 


f y [om] 
0 S 10 15 20 25 30 
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Fig. 14. Comparison of induced velocities obtained from the tests 
on the models with the original stern and modified one . 


Therefore the experimental tests fully confirmed 
the put hypothesis as follows : 


Calculations of the velocities around ship hull, induced by 
whirl systems representing the propeller itself and propeller 
race, make it necessary to modify input data to Biot-Savart 
formula, and as a result of the tests in question the proposed 

modification of the input data has been proved correct. 


TESTS ON FULL-SCALE SHIP 


The scientific aim of the presented investigations is to 
improve the algorithm applied in the software for calculating 
3D velocity field in the stern part of full-scale ship with taking 
into account propeller operation [7]. Therefore an important 
element of the investigations is to verify such field on a full- 
-scale ship. 

It is very hard to achieve reliable results from full-scale 
measurements of such field. In the subject-matter literature are 
known results of the measurements performed, both in model 
- and full-scale, on the hull of HSVA tanker ship, realized under 
the auspices of the model basin in Hamburg. 

The measurements of the velocity field before the propel- 
ler working on the full-scale ship were carried out through 
a window panel fitted in the stern part of the hull. They were 
conducted with the use of a laser anemometer but only respec- 
tive to the axial component of total velocity (i.e. with taking 
into account propeller-induced velocities) The measurement 


results are presented in the form of diagrams of isotachs, and 
only in the range covered by laser beam (Fig.15). The measu- 
rements were performed at the distance X = 0.21D from aft 
perpendicular. The broken line denotes the propeller circle of 
the diameter D = 6.1 m. 


(Hamburg Test Case) : full scale total axial 
velocity distribution; Vs = 9.26 m/s 
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Fig.15. Results of the measurements of axial velocity component, 
performed on the ship with operating propeller . 


And, in Fig.16 are presented the calculation results obta- 
ined from the modified PANSHIP software, also concerning 
only the axial component of velocity in the same cross-section 
before the propeller. 
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Axial component of effective velocity field behind the hull, V,/V, 
Effective wake fraction w, = 0.2249[-] 
Effective wake fraction w, , = 0.2273[-] 


Fig. 16. Results of the measurements of axial velocity component, 
performed on the ship with operating propeller . 


On the basis of analysis of the achieved results of calcu- 
lations, performed on the background of measurement results 
(compare Fig.15 and Fig.16), a qualitative similarity of both 
fields can be stated. Quantitative comparison can be more 
clearly presented in another form. In Fig.17, 18 and 19 the 
same results are shown in the form of diagrams of velocity at 
a given radius. To this end three radiuses : r/R = 1.0, 0.7 and 
0.5 were selected. 


If only accuracy of the measurements on the full-scale 
ship are taken into consideration (the curves presented in the 
figures should be symmetrical respective to the ship plane of 
symmetry) then the so-presented results are found astonishin- 
gly similar for the radiuses r/R = 1.0, and especially r/R = 0.7). 
It means that the PANSHIP software correctly determines the 
velocity field in the behind- the- hull flow and correctly expres- 
ses the velocity field induced by the propeller. In the diagrams 
are presented total velocity values which are formed along 
a considerable length of hull stern part at a significant share 
of propeller-induced velocities. It confirms that the PANSHIP 
can be successfully applied to the scaling of velocity fields on 
full-scale ship [7]. 
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Fig. 17. Comparison of the axial component of before-the-propeller 
velocity at the radius r/R = 1.0, obtained from measurements 
and calculations for the full-scale ship, respectively . 
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Fig. 18. Comparison of the axial component of before- the - propeller 
velocity at the radius r/R = 0.7, obtained from measurements 
and calculations for the full-scale ship, respectively . 


Calculation 
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Fig. 19. Comparison of the axial component of before- the - propeller 
velocity at the radius r/R = 0.5, obtained from measurements 
and calculations for the full-scale ship, respectively . 


FINAL REMARKS 


O It can be concluded that the obtained experimental tests 
fully confirmed the proposed hypothesis : 
Calculations of velocities around the ship, induced by 
whirl systems representing the propeller itself and pro- 
peller race require input data to Biot-Savart formula 
to be modified. 


-50 0 
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O The proved hypothesis can be considered as a kind of disco- 
very in fluid mechanics (in handbooks on hydromechanics 
and subject-matter literature no mention on that theme can 
be found) 


O However further theoretical research aimed at building 
a correct form of Biot-Savart formula in multiply connected 
space, is necessary 


O The proposed modification of the way of calculations 
of velocities induced around ship hull can be tentatively 
implemented in engineering practice as an effective appro- 
ximation. 


NOMENCLATURE 


— hull breadth 
— propeller diameter 
— hull length 
— number of propeller revolutions per second 
— radius of cylindrical cross-section around propeller axis 
— propeller radius 
> — hull resistance 
— resistance of hull with operating propeller 
— thrust deduction 
— propeller thrust, also hull draught 
— model speed 
ship speed 


: 2 2 2 
— total velocity y (Vut W) +V, +V, 


War rgow 


Ei 


ns 


(e) 


— induced velocity 


— total induced velocity 4/ v? + vy + v? 


x, Vy, Vz — velocity components 
Xp eZ — coordinates of propeller axis location 
y — distance of measurement points from hull plane of symmetry 


ô — hull block coefficient 

A — hull buoyancy 

AV,— velocity induced by whirl filament’s element 
à — model scale 


ro) 


<< < < eH 
| 
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a C onference —_ 
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of the Section 
on Exploitation 
Foundations 


On 25 May 2006 the Regional Group of the Section on 
Exploitation Foundations, Machine Building Committee, 
Polish Academy of Sciences (PAS), held its successive 
scientific seminar organized by Faculty of Engineering 
Sciences , Warmia - Mazury University in Olsztyn. 


Scientific workers of the Faculty 
presented the following papers : 


= A method for improving operation processes of track 
engine — by B. Kolator 

= Application of Exsys Covrid to maintenance of machi- 
nes — by K. Ligier and A. Rychlik 

= A system for maintaining the machines in tero-techno- 
logical approach — by P. Mikołajczak 


After discussion and replies from the side of the authors 
to questions directed to them, the organizers presented 
scientific laboratories of the Faculty. 


Miscellanea 


The books Ship Turbine Power Plants: Fundamentals of Thermodynamical Cycles and Introduction to the Theory of 
Marine Turbines by Krzysztof Kosowski form the first two volumes of a series of four monographs on marine turbine 
power plants. They are meant for mechanical engineers and for graduate students of technical universities, as well as 
marine and naval academies. They were initially elaborated as part of the EuroMTEC program for the module Advanced 
Ship Propulsion and Equipment, and have now been remarkably developed and noticeably extended. When writing them 
the author made use of some of the most outstanding works on the subject and to take into account his experiences from 
work at university. Fundamentals of power plant cycles and principles of turbine operation were laid out using renowned 
books ranging from the first works of A. Stodola published 100 years ago to the latest scientific papers and information 
given by major turbine producers. 


Ship Turbine Power Plants: Fundamentals of Thermodynamical Cycles 
This book deals with thermodynamical cycles of steam and gas turbines, and turbine power plant 
arrangements. 


Chapter 1: Fundamental Principles of Thermodynamics 
Chapter 2: Steam Turbine Cycles 

Chapter 3: Gas Turbine Cycles 

Chapter 4: Combined Turbine Cycles 


(ISBN: 83-922007-2-1, Published by: Foundation for the Promotion of Maritime Industry, Gdansk, 
2005, 280 pages, hardback, full colour) 


Introduction to the Theory of Marine Turbines 
This book deals with the fundamental aspects of axial turbine theory. 


Chapter 1: Review of Gas Dynamics 

Chapter 2: Axial Stage Theory 

Chapter 3: Stage Internal Losses 

Chapter 4: Efficiency Characteristics 

Chapter 5: Calculations of Flow in Turbine Stages 
Chapter 6: Multi-stage Turbines 


(ISBN: 83-922007-3-X, Published by: Foundation for the Promotion of Maritime Industry, Gdansk, 
2005, 261 pages, hardback, full colour) 


The books have earned high appreciation of the reviewers: 

„The author deals not only with turbine theory, but also with a wide range of aspects connected with strength, opera- 
tion, technology and dynamics. [...] The discussed problems are presented in a clear and concise way [...] and the book 
meets all standards of an academic handbook.” 


„The handbook contains not only a classical approach to the principles of turbomachinery, but also presents state-of- 
-the-art design methods and results of recent research.” 


„The form of presentation of the material and the collected examples deserve admiration [...], while the text and 
figures are excellently chosen.” 
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Analysis of hull resistance 
of pushed barges in shallow water 
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Jan Kulczyk 

Maciej Zawislak 

Wroctaw University of Technology 


ABSTRACT 


These authors performed a set of numerical calculations of water flow around pushed barges differing to 

each other by bow forms. The calculations were executed by means of FLUENT computer software. Tur- 

bulent free-surface flow of viscous liquid was considered. In this paper the calculated values of barge hull 
resistance split into bow, cylindrical and stern part components, have been compared and presented. 


Keywords : inland waterways ship, hull resistance 


INTRODUCTION 


The barges operating in push-train mode are characterized 
by great values of the hull block coefficients (C> 0.85), that 
ensures achieving large values of their displacement at assumed 
main dimensions. On the other hand, to decrease their building 
costs, usually is applied a simplified bow form consisted of 
practically developable surfaces divided by chine lines, thus 
relatively simple in building. Such approach is a rational and 
economical compromise since service speed of ships on inland 
waterways is of the order of 10 — 15 km/h. 

An inspiration to undertake the research on hull resistance 
of inland navigation cargo ships has been given by the informa- 
tion coming from an inland navigation ship owner that the fuel 
consumption on a given shipping route differs significantly in 
the case of pushed barges differing to each other first of all by 
their bow forms. These authors decided to investigate which 
bow forms of pushed barges ensure obtaining the smallest hull 
resistance values. To this end several characteristic bow forms 
were selected [1]. Each of the selected characteristic forms 
has been adjusted to a barge having the main dimensions : 
Lo x B x T= 48.75 x 9.0 x 1.7 m, under the assumption that 
the bow length from the end of the cylindrical midship body up 
to the bow transom plane is equal to L, = 8.0 m. Next, series 
of calculations of the flow around the push-trains consisted of 
two barges connected to each other by their stern parts, were 
performed. The calculations were executed by means of the 
FLUENT commercial computer software which makes it pos- 
sible to take into consideration all factors of crucial influence 
on ship resistance, i.e. viscosity of water, turbulence of flow, as 
well as wave system on water free-surface around the ship. 

Quality of the calculation results of free-surface water flow 
around inland navigation ships in shallow water, has been as- 
sessed during the previous research investigations carried out 
by these authors [2 , 3]. In view of a limited performance of 
the computers being at the authors’ disposal most of the com- 
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putations was performed for the hulls in a reduced scale. As 
a rule the same scale has been applied as in the case of model 
testing in a towing tank. A direct comparison of the results of 
the calculations with those from model tests has confirmed that 
the applied software is useful in calculating hull resistance and 
determining wave profile on the ship side. 


HULL FORMS 
OF THE CONSIDERED BARGES 


The calculations of water flow around hulls of the barges 
were performed for 11 trains of barges fitted with bows of the 
following forms (Fig.1) : EIIB, EBM, EUBV2, EIBH, ELI, 
ELIM, WALE, WALC, B, B3 and HEL. The first four constitute 
a group of similar forms. They have been elaborated on the basis 
of the hull form of the EUROPA IIb pushed barge popular on 
the West - European waterways. They differ to each other by 
the shape of longitudinal cross-section contour in the plane of 
symmetry. For the four barges similar results were achieved. 
The hull form of ELI barge has been elaborated on the basis of 
an elliptical bow (Ellipsenbug) proposed by Nussbaum [4]. The 
ELIM form is a simplified version of the ELI form. Rounded 
segments of frame sections have been replaced by straight-line 
segments inclined by the angle of 45°. As a result, the surface 
between the bottom and side of hull has become a developable 
surface. The WALE and WALC forms have vertical sides and 
are of the simplest geometry. They differ to each other only by 
a shape of water-planes which are elliptical in the first case, 
and in the other — circular segments tangent to ship sides. 
The B and B3 bow forms have been designed by the team 
working on the project. The form B ensures obtaining a high 
block coefficient value of the bow. Owing to the flat bottom 
it is possible to make the barge cubicoid hold much longer. 
The other bow form is more fine — it has a more inclined stem 
and higher elevated line of the side chine. The HEL bow form 
has been designed by these authors. Side surface of the bow 
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PP S PP 
EBM | | 
KLW 
2 2 
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KLW ] 
g er 
pp |} | o 
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HEL B3 
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Fig. 1. Hull forms of the considered barges. PS - Plane of symmetry, PP - Base Plane, KLW - Design waterline . 


is a fragment of a regular helicoid with its axis located in the 
plane of hull side. This shape was assumed to jostle water aside 
like the WALC form and simultaneously to integrate a broad 
deck and bow transom like in the case of the EIB bow form. 
In the below presented table, are given values of the block 
coefficient of hull and that of bow which has been defined as 
follows : Cpe = V,/(L, x B x T). 
The authors have intended to check if any unambiguous 
relation between values of the above mentioned coefficients 
and hull resistance, takes place. 


EIB 
ENUBM 
ENBV2 
EIBH 
ELI 


ELIM 

WALE 

WALC 
B 


B3 0.952 0.663 
HEL 0.939 0.589 


The calculations were performed for the train of two barges 
connected to each other by their stern parts. To elaborate grids 
for numerical calculations the assumption was made that the 
stern form of a single barge influences train’s resistance to the 
same degree, irrespective of an applied bow version. For this 
reason the aft bottom undercut was not modelled and the 
cylindrical parts of both barges were made longer and joined 
together in the aft transom plane. 

The identical flat bilge form of 200 mm in height, was ap- 
plied to all the barges, except of those having ELI bow form, 
where the cylindrical bilge form of 200 mm radius was used. 


CALCULATION CONDITIONS 


The flow calculations were performed for two values of 
water depth: 2.0 m and 3.4 m. The first of them models the 
conditions of very shallow water (h/T = 1.18). In the case of 
canalised rivers such conditions appear only in certain places 
— along short sections of a waterway. The other water depth 
(h = 3.4m) models shallow water conditions, which is more real- 
istic for average service conditions on the domestic waterways. 

In both the cases the calculations were carried out for the 
ship speed equal to 2.48 m/s (8.93 km/h). i.e. at the Froude 
number Fn, = 0.56 in a more shallower water, and Fn, = 0.43 
in a deeper water. At such speed a significant sagging of the 
ship should be taken into account, especially at the water depth 
equal to 2.0 m. However the taking of sagging into account in 
calculations is associated with a change of location of bound- 
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aries of computation area and a significant increase of time of 
computations. The authors have assumed that the neglecting of 
sagging introduces the same errors to resistance values of all 
the considered hull forms. Hence the differences in calculated 
resistance values would maintain the same, and to compare 
directly the bow forms would be possible. 

For the calculations performed within the frame of this 
research work the authors made use of the same principles of 
building the computational grids and controlling calculation 
runs as those used in the previous research work [3]. 

All the calculations were performed in the model-scale 
of 1:14. The computational grid covered the rectangular area 
extending up to 41.25 m ahead the bow and behind the stern, 
and 45.5 m overboard. The grid mesh was so designed as to 
ensure precise modelling the hull form and the flow around hull 
surface. As a rule a regular grid consisted of cubicoid elements 
was applied, but irregular one — only locally. For the reason of 
a limited computer performance the number of elements did 
not exceed 200 000. 

In the FLUENT software the problem in question was defined 
as non-stationary one. The equations were integrated till reach- 
ing a stationary state. The applied time-step of 0.01s ensured 
reaching the convergence of calculations after 30 000 steps. 
To model the turbulence phenomenon the model RNG k-e was 
selected. A single run of calculations took 48 h on average. 


RESULTS 


In contrast to the experimental methods (towing tank model 
tests) the application of the numerical computation method to 
fluid mechanics (CFD) makes it possible to split hull resistance 
into the components resulting from normal stresses (pressu- 
re) and tangential stresses (liquid viscosity). It also makes it 


Total resistance 


possible to easily determine forces acting on various parts of 
the hull. In order to analyse a contribution of particular hull 
segments in total hull resistance the authors split the entire hull 
surface into three parts : 


** the bow (from the bow transom plane of fore barge to the 
cylindrical midship body) 

% the cylindrical midship body (precisely — joined midship 
bodies of both fore and aft barges), and 

** the stern (from the cylindrical midship body of aft barge to 

the bow transom plane of aft barge) (Fig.2). 


Vs 
—- 


81.5 m 


8.0m 
T 1 
Cylindrical midship body 
Bow 


Fig. 2. The train of two barges split into three parts 
for resistance analysis purposes 


The below presented values of hull resistance, calculated 
by means of the FLUENT software for free-surface flow con- 
ditions, take into account hydrostatic pressure. 

In further considerations it was assumed that hull resistance 
is a force acting in the direction opposite to ship speed vector 
(i.e. aft). A negative value of stern resistance means that the 
resultant force acting onto the stern is directed fore. The greater 
the force the smaller the total hull resistance. In Fig.3 and 4 the 
bow forms are ranked in a sequence resulting from increasing 
value of total resistance. 

In design and service practice, quality of a pushed barge hull 
form is assessed by using the unit resistance values, 1. e. those 
taken per unit buoyancy or volume of underwater part of ship’s 
hull. The unit resistance values are compared in Fig. 5. 


Bow resistance 
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Fig. 3. Resistance of the two-barge-train model at the water depth of 2.0m . 
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Fig. 4. Resistance of the two-barge-train model at the water depth of 3.4 m . 
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is exerted a greater pressure resistance and smaller viscosity 


resistance than onto the remaining hulls (Fig. 3 and 4). Fig. 6. The bow block coefficient C, (the bows are ranked on the basis 


of their hull resistance values at the water depth of 3.4 m, see Fig.4) . 
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© In Fig. 7 and 8 the points are clustered in two groups. The bow forms : HEL, B, WALC and WALE belong to the first group, 
the remaining — to the other group. The bows of the first group have a straight, vertical or only slightly inclined stem, and 
greater resistance values as well. This observation suggests that the vertical or only slightly inclined stem is not favourable 
from the point of view of pushed barge hull resistance. 
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Fig. 7. The relationship of hull resistance and the hull block coefficient C, . 
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Fig. 8. The relationship of hull resistance and the bow block coefficient C,,, . 


NOMENCLATURE 
B - hull breadth 
C,  - hull block coefficient 
C - bow block coefficient 
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RVM - viscosity resistance in model-scale 
RMM- resistance of cylindrical midship body in model-scale 
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onference 
—Corterence 
HYDROACOUSTICS 2006 


On 23 — 26 May 2006 
at Krynica Morska upon Vistula Bay was held : 


13 SYMPOSIUM ON HYDROACOUSTICS 


organized by the Department of Marine Electronic Sys- 

tems, Faculty of Electronics, Telecommunication and 

Informatics, Gdansk University of Technology, under the 

auspices of : European Acoustics Association, Hydro- 

-acoustics Section of Committee on Acoustics, Polish 

Academy of Sciences, and Gdansk Division of Polish 
Acoustical Society. 


The Symposium was commenced 
by the key-note lecture on : 


Research and development on underwater acoustic sys- 
tems of Polish Naval University and Gdansk University of 
Technology for the Polish Navy — by G. Grelowska (Polish 
Naval University) and L. Kilian (Gdansk University of 
Technology). 


During 4 plenary session of the Symposium 
the following 5 invited papers were presented : 


> Science and technology in Polish Ministry of Defense 
by W. Drag (Polish Ministry of Defense) 

> New scientific multi-beam systems for fishery research 
applications — by L. Nonboe (SIMRAD, Norway) 

> The state of the Baltic Sea hydro-acoustical investiga- 
tions (selected problems) — by Z. Klusek (Institute of 
Oceanology, Polish Academy of Sciences) 

> Synthesis and wavelet analysis of side-scan sonar sea 
bottom imagery — by J. Tegowski (Institute of Oceano- 
logy, Polish Academy of Sciences) and A. Zielinski 
(University of Victoria, Canada) 

> Quadrature phase detection in an acoustic positio- 
ning system — by A. Zielinski (University of Victoria, 
Canada) and Y.Shi (Southwest Jiaotong University, 
China) 


The remaining 25 papers were presented during 4 
panel sessions. Original papers, both theoretical and 
experimental, concerning problems of hydro-acoustics 

and its applications are published in the annual journal 
Hydro-acoustics. e 


Miscellanea 
Workshops 2006 


Under this name, on 30 March + 1 April 2006, Faculty 
of Maritime Technology, Technical University of Szcze- 
cin, arranged the series of popular scientific lectures and 

demonstrations to promote the courses on 


Ocean Engineering and Transport 
conducted at the Faculty. 
Academic lecturers presented the following themes : 


* Safety at sea — by M. Hann 

* Gas an oil mining from sea bed 
by W. Chadzynski 

* Contemporary maritime industry 
and shipping — by T. Jastrzebski 

* Neural networks — by D. Pielka 

* Shapes of sound — by S. Weyna 

* =Super-computers and turbulence 
by T. Abramowski 

* Unconventional energy sources on ships 
by W. Zenczak 

* Digital evolution — by P. Nikonczuk 

* Stirling’ engine — by A. Żmuda 


The last day the underwater apparatuses built at the 
Faculty were demonstrated. The Workshops appeared 
very interesting for many visitors hence it was decided to 

organize them every year. 


oe J) 


f+ C onference —* 


SEM — ECO 


On 12 May 2006 the scientific seminar on : Ecological 
problems in operation of combustion engines, organized by 
Prof. L. Piaseczny, was held at Polish Naval University. 

The seminar program contained two lectures presented 
by the scientific workers from Warsaw University of 
Technology, namely : 


# Selected problems of emission of PM10 solid partic- 
les from exhaust gas systems of combustion engines 
by M. Zegota 

#* Development trends of combustion engines for usage 
vehicles — by Z. Chtopek 


Both the topics triggered very interesting discussion 
which enriched the knowledge passed on in the presented 
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Stiffness characteristics and thermal 
deformations of the frame 
of high power marine engine 


Lech Murawski 

Marek Szmyt 

Centrum Techniki Okrętowej S. A. (CTO) 
(Ship Design and Research Centre) 


ABSTRACT 


In the subject-matter literature detail data on stiffness of the crankshaft foundation connected with the 
frame of marine main engine are still lacking. Thermal deformation models of the engine 5 casing, propo- 
sed by engine producers, are excessively simplified. However the parameters are crucial for the shaft-line 
alignment analysis as well as for the analysis of interactions between the shaft-line and engine crankshaft, 
especially in the case of high power engines. This paper presents a determination method of the marine 
engine body characteristics as well as results of example computations performed for a Sulzer 7 RTA 84 C 
engine installed on a~3000 TEU container ship. It has been demonstrated that the producers assumption 
about parallel displacement of the crankshaft axis in thermal working conditions is too rough. The ther- 
mal deformation of the engine is of hogging character, which results in significant change of the moment 
load exerted on the crankshaft and shaft line. The stiffness parameters recommended by the producers 
for the shaft-line alignment are estimated correctly, however they represent only engine body flexibility, 
without taking into account ship's hull flexibility. 


Keywords : marine main engine, main bearing, static and dynamic stiffness characteristics of bearing, 
thermal deformation, temperature distribution 


INTRODUCTION 


The aim of the presented work has been to evaluate dis- 
placements of the engine crankshaft axis in different working 
conditions of propulsion system [3 , 5]. In the shaft-line ali- 
gnment methods the crankshaft-shaft line interaction has been 
considered in a simplified way so far. The crankshaft has been 
modelled as a linear system of cylindrical beam elements, 
whereas its thermal displacements and its foundation stiffness 
have been evaluated on the basis of simple data supplied by the 
producer without considering the type of the ship on which the 
engine has to be installed [4]. The goal of the work has been 
to improve representation of the boundary conditions of the 
marine power transmission system. It is especially important 
for the high power propulsion systems since in the literature 
many examples of the damage of the first three main bearings 
(counting from the driving end) of the main engine can be 
found [1 , 2]. One of the possible causes of such state might 
be the imprecise mathematical model of crankshaft, proposed 
for the analysis of shaft line alignment. 

Within the frame of this work several analyses of the Sul- 
zer 7 RTA 84C engine installed on a big container ship (of 
~3000 TEU capacity) were carried out. Also, the computation 
of engine body deformation under gravity load as well as 
the analysis of its thermal deformation in nominal working 
conditions was performed. The static stiffness (horizontal and 
vertical) of each of the main bearings were evaluated and then 
their dynamic stiffness was determined in the frequency range 
of 0+30 Hz. As the forced vibration analysis was performed 
with the use of the modal superposition method, it was ne- 
cessary to determine in advance the natural frequencies and 
eigenvalues in the frequency range of 0+70 Hz. The thermal 
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analysis requires an accurate temperature distribution on the 
engine body to be known. Appropriate data were obtained from 
comprehensive temperature measurements performed on the 
ship and her main engine. 


ANALYSIS METHOD 


The FEM model of the body of Sulzer 7 RTA 84 C engine 
is presented in Fig. 1. Fig.2 shows a part of the model represen- 
ting the engine main bearing. The engine body model contains 
almost 200 thousand plate and solid elements of over 930 
thousand degrees of freedom. 

The subject of the analysis was the body of Sulzer 
7 RTA 84 C engine. The analysis involved, apart from heat flow, 
also thermal deformation and stress calculated by means of a 3D 
heat transfer model. The method is based on the solution of the 
heat flow equation (with variable coefficients). For stationary 
heat flow the equation has the following form: 


& k(x. y.2) + o kosy, 4 D 
+Ê kx, y,z)3L = (0) 


If the convective boundary conditions on both hot and cold 
surface are assumed then the following relations are valid : 


-k(x,y,z)Z = h(T,-T) (2a) 


-k(x, y,z)ĝ = h(T,-T,) (2b) 


y 
Pus 


Fig. 1. FEM model of Sulzer 7 RTA 84 C engine . 


y 
Aa 


Fig. 2. Main bearing frame of Sulzer 7 RTA 84 C engine . 


No heat flow is assumed on other surfaces of boundary 
planes of the considered hull segment, which leads to the 
equation: 

oT _ 
ay (3) 

The FEM thermal analysis was performed with the use of 
MSC NASTRAN software. MSC PATRAN software was used 
as a pre-and post-processor for the calculation results of the 
stress under thermal load. 

Before starting the thermal deformation analysis of the 
engine it is necessary to determine temperature distribution on 
its body. The temperature map was created on the basis of the 
measurements carried out on a marine main whose size and 
structure was similar to those of the engine in question. 


TEMPERATURE DISTRIBUTION 
MEASUREMENT ON ENGINE BODY 


The temperature distribution measurements on the engine 
body were performed on the Sulzer 8 RTA 96 C engine during 
sea trials. The engine load was kept stable in nominal working 
conditions. Alfa-Tech Rytek MT 4 pyrometer was used for the 
measurements. The example layout of measurement points (on 
the port side) is shown in Fig.3. The results of measurements 
in those points are presented in Tab. 1. On the starboard side as 
well as on the fore and aft end of the engine the measurement 
points were distributed in a similar way. 


*13 


"S 179 


Fig. 3. Layout of the measurement points on the port side 
of the main engine body . 


Tab. 1. Results of the temperature measurements performed 
on the port side of the main engine body . 


Port side of the Engine 


Measured 
tempera- 
ture [°C] 


Measured | Measure- 
tempera- ment 
ture [°C] | point No. 


Measured | Measure- 
tempera- ment 
ture [°C] | point No. 


Measure- 
ment 
point No. 


THERMAL ANALYSIS 
OF ENGINE BODY DEFORMATION 


In the numerical thermal analysis the value of heat conduc- 
tivity coefficient for steel was assumed equal to 42.9 [W/m K]. 
The heat flow analysis was performed for the condition of the 
hold being hot and thermally balanced. Because of lack of more 
precise data the heat transfer coefficients were assumed as for 
heated cargo in accordance with DNV Classification Rules. 
The values of heat transfer coefficients are presented in Tab.2. 
The thermal expansion coefficient of the engine body was 
assumed equal to a = 1.6x10°. The temperature distribution 
on the engine body, shown (in °C) in Fig.4, was analogous to 
that obtained from the measurements. 


Tab. 2. Assumed values of heat transfer coefficients [W/m °C] . 


From air in the hold to the inner bottom structure 58.1 
From air in the hold to the side structure 58.1 
From air in the hold to the deck structure 58.1 

From outboard water to the hull shell 7400 
From air to the hull shell 23.2 
From air in the double bottom to the hull structure 0 
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y 


ar 


Fig. 4. Temperature distribution on the engine body, assumed 
for the analysis . Note : On this and all next figures 
the SI standard units (e.g. m, Pa) are applied . 


In Fig.5 the thermal deformation of the engine body is 
presented in the form of fringe plot. From the point of view of 
the propulsion system and main engine — shaft line interaction, 
displacements of the engine main bearings are most important. 
The values of the displacements are presented in Tab.3. The 
bearings are numbered beginning from the driving end of the 
shaft line (the right-hand side of Fig.5). 


mx 


Fig. 5. Thermal deformation of the body of Sulzer 7 RTA 84 C engine . 
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Tab.3. Thermal displacements of the main bearings 
of Sulzer 7 RTA 84 C engine . 


Vertical 
displacement 
[mm] 
0.756 


Axial 
displacement 
[mm] 
1.905 


Horizontal 
displacement 
[mm] 
-0.004 


Main bearing 
No. 


2 
3 
4 
5 
6 
7 
8 
9 


For the crankshaft axis translation the main engine’s 
producer recommends to use the following formula : 


h. = dh(T, - Ta) (4) 
where : 
h, — displacement of the main bearing axis [mm] 
dh — thermal expansion coefficient of the engine 
body [mm/K] 
T, — service temperature of the engine [K] 
T, — air temperature in the engine room [K]. 
For the examined Sulzer 7 RTA 84 C engine 
the translation of the main bearing axis amounts to : 
h,=18.4-107°(53.0-20)=0.607mm (5) 


The numerically computed value of the translation of the 
shaft line axis (Tab. 3) is greater than that recommended by 
the producer despite the fact that the measured temperature 
of the real engine is slightly lower than its specified service 
temperature. The difference for the first bearing is not par- 
ticularly large (less than 20%), but other bearings are much 
more displaced. It seems that the producer’s assumption on 
the parallel translation of the crankshaft axis is incorrect. The 
hogging deformation of the crankshaft results in a significant 
change of the moment transferred from the shaft line. Its effect 
calculated by using a precise shaft line alignment analysis 
would be considerable. 


ANALYSIS OF STATIC STIFFNESS 
OF MAIN BEARINGS OF THE ENGINE 


Determination of the static stiffness consists in applying 
unitary forces equal to mass forces and radial gas forces (of 
750 KN), to each ofthe main bearings, one by one, first in verti- 
cal and then in horizontal direction. The achieved displacements 
serve to calculate the local static stiffness. The quantities are 
very important for shaft line alignment analysis since applying 
only the ship’s hull stiffness may be insufficient. The main 
engine producers usually provide (on request) information on 
the crankshaft foundation stiffness but without dividing the 
quantity into that concerning engine body alone and ship hull. 
Some time ago this parameter was assumed infinitely large, 
now it is considered as large as 6.0x10° N/m. The CTO Co. 
gained vast experience concerning the stiffness of ship hulls 
of many types but no main engine body stiffness has been so 
far examined by this company. 

The deformation of the main engine body under only gravity 
load was computed first. It was observed that the deformation 
due to gravity load was several times smaller (of the order of 
only 0.01mm) than that due to thermal load, cylinder mass and 


gas forces. Therefore in further analysis the influence of gravity 
load may be neglected. 

Next, was carried out the static analysis consisting in 
application of the load first in vertical and then in horizontal 
direction (18 load cases). The estimated radial forces were 
distributed on the main bearing surface. Tab.4 contains the 
static stiffness values for particular main bearings (numbered 
beginning from the aft — driving end of the crankshaft). The 
selected deformation values and stress fringe plots are pre- 
sented in Fig. 6+8. 


Tab. 4. Static stiffness of the main bearings of Sulzer 7 RTA 84 C engine . 


Main bearing No. 


Horizontal stiffness 
x10° [N/m] 


10.942 
11.232 


Vertical stiffness 
x10? [N/m] 


6.615 
6.828 


He 


Fig. 6. Deformation of the main engine body 
under vertical load applied to 6 main bearing . 


Fig. 7. Deformation of the main engine body 
under horizontal load applied to 1“ main bearing . 


Fig. 8. Stress in the main engine body under vertical 
and horizontal load applied to 3" main bearing . 


On the basis of the performed analysis it can be stated 
that the static stiffness values specified by the producer are 
properly evaluated, however they represent only the flexibility 
of the engine body without taking into account the ship hull 
flexibility. In the authors’ opinion in the analysis of shaft line 
alignment the sum of both flexibility parameters should be taken 
into account. It can be observed that the engine body stiffness 
is very high and the cylinder mass and gas forces acting on 
a single cylinder have a little influence on displacements of 
other main bearings. 

It means that it is not required to take into account any 
coupling between particular bearings — hence there is no neces- 
sity to determine the equivalent stiffness reflecting the engine 
structure integrity. It is intriguing that the horizontal stiffnesses 
are higher than the vertical ones, not so as in the case of most 
other marine structures. The stress level in the main bearing 
structure is not high as it does not exceed 15 MPa for horizontal 
load and 22 MPa for vertical one. 


ANALYSIS OF DYNAMIC STIFFNESS 
CHARACTERISTICS OF BEARINGS 


The determination of the dynamic stiffness characteristics 
consists in applying unitary forces defined as a function of 
excitation frequency, to each of the main bearings, one by one, 
first in vertical and then in horizontal direction. The obtained 
displacements serve to calculate the local dynamic stiffness for 


POLISH MARITIME RESEARCH, No 1/2007 19 


OPERATION & ECONOMY 


OPERATION & ECONOMY 


each excitation frequency, independently. These quantities are 
important for shaft line lateral vibration analysis. The data on 
the dynamic stiffness of main engine body are not available 
neither from literature nor from the producers. 

The analysis of the dynamic stiffness of Sulzer 7 RTA 84 C 
engine body was carried out by using the MSC Nastran FEM 
solver and modal superposition method. The stiffness characte- 
ristics were determined in the frequency range of 0+30 Hz; this 
way the full spectrum of possible propulsion system excitation 
frequencies was covered. In such case according to the authors' 
experience, it is recommended to calculate normal modes in 
the range of natural frequency, taken at least twice as wide as 
that. The eigenmodes were determined for frequency values up 
to 70 Hz. The eigenvector maps for the most important normal 
modes are shown in Fig. 9+11. The applied denotation for engi- 
ne body deformations, commonly used by engine producers, is 
as follows: H — lateral deformation, X — twisting deformation, 
C — bending deformation in vertical plane. 


Fig. 10. ”X”- type normal mode at the natural frequency of 45.30 Hz . 


It is important that the first significant normal modes have 
their natural frequencies above the range of excitation frequ- 
encies of the propulsion system. And, the significant normal 
modes are only a few and those of interest concern the whole 
engine body. No significant normal modes were found in the 
region of the engine main bearings. It speaks well for the correct 
design of the engine body, i.e. of sufficiently rigid structure. 
In such case the resonance discontinuities of main bearing 


20 POLISH MARITIME RESEARCH, No 1/2007 


J 


Fig. 11. ”C”- type normal mode at the natural frequency of 69.95 Hz . 


flexibility are not expected and the characteristics appear close 
to linear. For this type of structure it is allowed to use only 
static stiffness as the expected dynamic amplification should 
be insignificant. 

Next step was to compute forced vibration by applying first 
vertical load and then horizontal one to each of the engine main 
bearings, one by one (18 load cases). The analysis was performed 
within the frequency range of 0+30 Hz. Tab.5 and 6 contain 
dynamic stiffness values for characteristic excitation frequencies 
and particular bearings (numbered beginning from the aft — driv- 
ing end of the crankshaft). The nominal rotational speed of the 
examined engine was 100 rpm. The engine had seven cylinders, 
and the propeller - five blades. For such configuration the basic 
excitation frequencies were 8.33 Hz and 11.67 Hz. 


Tab. 5. Vertical stiffness of main bearings of Sulzer 7 RTA 84 C engine . 


Stiffness 
values 
at 0 Hz 

x10° [N/m] 


Stiffness 
values 
at 8.33 Hz 
x10° [N/m] 


Stiffness 
values values 
at 11.67 Hz| at 25 Hz 
10° [N/m] | x10° [N/m] 


Main Stiffness 


Bearing 
No. 


niana 


Tab. 6. Horizontal stiffness of main bearings of Sulzer 7 RTA 84 C engine . 


Stiffness 
values 
at 0 Hz 

x10° [N/m] 


Stiffness 
values 
at 8.33 Hz 
x10° [N/m] 


Stiffness 
values values 
at 11.67 Hz| at 25 Hz 
10° [N/m] | x10° [N/m] 


10.878 10.447 
10.616 


Main Stiffness 


bearing 


11.155 


The sample plots (for the main bearing No.3) of displacement of main bearing foundation in function of vertical (and subsequ- 
ently horizontal) frequency are presented in Fig.12. Selected fringe plots of the engine body deformation in the region of the 
main bearing No. 3 for excitations in both directions are shown in Fig.13. 


Displacement translational 


T T 1 
0 5.00+00 1.00+01 1.50+01 2.00+01 2.50+01 3.00+01 
Frequency 


0 
-1.50-05 
-3.00-05 
-4.50-05 
-6.00-05 


pe 


-7.50-05 


Displacement translational 


-9.00-054 1 i i 
0 5.00+00 1.00+01 1.50+01 2.00+01 2.50+01 3.00+01 
Frequency 


Fig. 12. Displacement of the main bearing No. 3 in function of excitation frequency . 


Fig. 13. Deformation of the main engine body under vertical and horizontal excitation applied to the main bearing No. 3. at the frequency of 8.33 Hz. 


The dynamic stiffness values for the excitation frequency 
of 0 Hz are almost identical with the static stiffness (their 
differences are at 3 decimal place — see Tab. 4, 5 and 6). It 
speaks well for the correctness of the dynamic analysis. As it 
was expected after performance the analysis of normal modes, 
the dynamic stiffness values did not significantly differ from 
the static ones. The stiffness decrease by 2% in vertical direc- 
tion may be observed and that in horizontal direction — by less 
than 7%. Such change of dynamic stiffnesses can not have any 
significant effect on the analysis of shaft line lateral vibration. 
In standard (commercial) analyses the stiffness evaluation 
may be limited to only a static quantity which can be assumed 
constant in the domain of excitation frequency. 


CONCLUSIONS 


O The numerically computed shaft line displacement is greater 
than that recommended by the producer despite the fact 
that the measured temperature of the real engine is slightly 
lower than its specified service temperature. For the first 
bearing (counting from the shaft line side) the difference 
is not particularly big (lower than 20%), but other bearings 


are much more displaced. Hence the producer’s assumption 
on the parallel translation of the crankshaft axis seems to 
be incorrect. The hogging deformation of the crankshaft 
results in a significant change of the moment load exerted 
by the shaft line. It may be expected that the effect result- 
ing from the precise shaft line alignment analysis would be 
considerable. 


O On the basis of performed analysis it can be stated that the 


static stiffness values specified by the producer are properly 
evaluated, however they represent only the engine body 
flexibility without taking into account the ship’s hull flex- 
ibility. In the shaft line alignment analysis the sum of both 
flexibility parameters should be taken into account. It can 
be observed that the engine’s body stiffness is very high 
and the cylinder mass and gas forces acting on one cylinder 
have a little influence on displacements of other main bear- 
ings. Therefore it is not necessary to take into account any 
coupling between particular bearings. 


O It is important that the first significant normal modes have 


natural frequencies above the range of excitation frequ- 
encies of the propulsion system. Moreover the significant 
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normal modes are few and those of interest concern the 
whole engine body. No significant normal modes were 
found in the region of the engine main bearings. It means 
that the engine body in question is a well-designed rigid 
structure. In such case the resonance discontinuities in the 
flexibility of main bearings should not be expected and their 
characteristics should be close to linear. 


O The dynamic stiffness values differ insignificantly from the 
static ones. The stiffness decrease by 2% can be observed 
in vertical direction and that by more than 7% in horizontal 
direction. Such change of dynamic stiffness values can not 
have any significant effect on results of the analysis of shaft 
line lateral vibration (whirling). In commercial analyses 
evaluation of stiffness may be limited only to determining 
the static value which — in the domain of excitation frequ- 
ency — can be assumed constant. 


O The direction of research in question looks very promising. 
It would make it possible to introduce such improvements 
to high power propulsion systems as to avoid failures of 
the engine's main bearings. The presented method may 
be also used for more advanced and complete numerical 
computations carried out for main engines of other types, 
installed on ships having specific hulls. 


O The next step in developing the proposed method of pro- 
pulsion system analysis should be incorporation of a more 
complex crankshaft representation based on its full 3D 
characteristics. The crankshaft springing effect on the shaft 
line alignment should be also examined. 
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NOMENCLATURE 


dh — thermal coefficient of expansion of the engine body [mm/K] 
h, — vertical displacement of the main bearing axis [mm] 

h —thermal heat transfer coefficient [W/(m°K)] 

k —thermal conduction coefficient [W/(mK)] 

T — temperature [K] 

T, — air temperature in the engine room [K] 

T, — service temperature of the engine [K]. 


BIBLIOGRAPHY 


1. MAN B&W Diesel A/S : Bearings. Copenhagen. 2000 

2. MAN B&W Diesel A/S : Elasto-hydro-dynamic evaluation of 
main bearing performance. Copenhagen. 2002 

3. American Bureau of Shipping: Guidance notes on propulsion 
shafting alignment. Houston. 2004 

4. MAN B&W Diesel A/S : Shafting alignment for direct coupled 
low-speed diesel propulsion plants. Copenhagen. 1995 

5. Wärtsilä: Sulzer RTA-C. Technology Review, Helsinki. 2003. 


CONTACT WITH THE AUTHORS 


Lech Murawski, D.Sc., M.E. 
Marek Szmyt, M.Sc., M.E. 
Centrum Techniki Okrętowej S. A. 
Rzeczypospolitej 8 
80-369 Gdańsk, POLAND 
e-mail : Lech.Murawski@cto.gda.pl 


22 POLISH MARITIME RESEARCH, No 1/2007 


Mourning 


9 August 2006 was the time of a deep sorrow 
for the circle of Polish shipbuilders as 


Professor Jerzy Doerffer 
passed away this day. 


Graduate of Shipbuilding Faculty of University 
of Glasgow and Gdańsk University of Technology. 
Professor of Gdańsk University of Technology 
where he worked since 1948. 

The organizer and the first Head of the Department 
of Ship Technology and its auxiliary unit. 

The Dean of Shipbuilding Faculty 
in the years 1953-54 and 1958-64. 

The Rector of Gdańsk University of Technology 
in the years 1964-67. 


The Chairman of the Forum of Shipbuilding 
and Ship Repair Industry. 


The creator of the scientific school in the domain 
of shipbuilding technology. 
He was a worldwide recognized authority 
on this domain. 


The co-author of novel design solutions 
and manufacturing techniques in ship technology. 


Doctor Honoris Causa of Gdańsk University of Tech- 

nology, Leningrad Shipbuilding Institute, University of 

Glasgow, University of Rostock, Polish Naval University 
and Technical University of Szczecin. 


He was also honoured with William Froude Medal 
by the Royal Institution of Naval Architects, London. 


Prof. Doerffer was not only an outstanding engineer-in- 

ventor but also a moral authority of wonderful personali- 

ty. Due to his attitude and work he won respect from the 
side of the whole circle of Polish shipbuilders. 


He was a figure-symbol of whom the circle 
of Polish shipbuilders boasts. 


Humanist and teacher. He promoted a whole genera- 
tion of Polish engineers in shipbuilding technology. 


A member of prestigious scientific 
institutions in Poland and in the world. 
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Usefulness assessment of standard measuring 
instruments installed on sea-going ships 
to perform energy measurements 


Zygmunt Gorski, 
Romuald Cwilewicz, 
Gdynia Maritime University 


ABSTRACT 


The presented work is a contribution to discussion on usefulness of application of measurement instrumen- 
tation used on sea-going ships for energy measurement and scientific research purposes. Contemporary 
sea-going ships are equipped as a rule with up-to-date measurement instrumentation usually based on 
electronic data processing and computer technique. These authors have made many times use of such 
instruments in their research work. This way it was not necessary to install any special instruments, that 
significantly reduced measurement cost. In such cases to obtain a sufficient accuracy of measurements 
constitutes a crucial problem. In this paper was presented an analysis of measurement errors of some 
operational parameters of ship and its main propulsion system, elaborated within the frame of the KBN 
research project no. 9 T12D 033 17. Results of the analysis confirm usefulness of the standard measurement 
instrumentation installed on ships, and its sufficient accuracy. 


Keywords : ship, standard measurement instrumentation, accuracy of measurements. 


INTRODUCTORY COMMENTS 


Contemporary sea-going ships are fitted as a rule witha rich 
set of control and measuring instruments. Commonly are in use 
analogous or digital measuring systems as well as those based 
on computer technique. They serve both for carrying out direct 
measurements and data processing in the systems for remote 
and programmed control of ship operation, main propulsion 
system and auxiliary shipboard devices. 

On ships a special high-accuracy measuring instrumentation 
is also applied. Such instrumentation is installed only for the 
time of measurements and it does not belong to ship’s standard 
equipment. 

Within the frame of the KBN research project no. 9 T12D 
033 17 some energy measurements were performed on the trai- 
ning-research ship Horyzont II in the period from 23.06.2001 
to 14.07.2001. The measurements covered the quantities 
associated with operation of ship main propulsion system in 
various service conditions, namely main engine (ME) power 
output, ship speed, fuel consumption by the main propulsion 
system. Additionally were measured selected parameters of ME 
operation, interesting from the point of view of loading state of 
main propulsion system, e.g. exhaust gas temperature, super- 
charging air pressure and temperature. Also, some parameters 
associated with ME operation such as exhaust gas content in 
various loading states, were investigated. For carrying out the 
measurements standard measuring instruments installed on 
the ship as well as special measuring devices installed by the 
research team, were used. 


The investigations were aimed among other at deter- 
mination of ship propulsion characteristics. To this end the 
measurements of ship speed, ME power output and its fuel 
consumption at various rotational speeds and propeller pitch 
settings, were performed. 

Specification of the standard measuring instruments installed 
on the ship and used for the research in question is as follows : 


> for the measuring of : propeller pitch setting, rotational 
speed of ME and its load index — the integrated computer- 
-based system for measuring operational parameters and 
control of settings of ship’s main propulsion system 
-Wichmatic 2 Propulsion Control System, Graphic System 
WM 11G ECR, WM2 12G ECR, WÄRTSILÄ NSD Norway 
AS, having digital read-out system 

> forthe measuring of ship speed : the SAL R1 U/N 701488C 
electro-magnetic log of Consilium Marine Co, having di- 
gital read-out system 

> for the measuring of : exhaust gas temperature at outlets 
from engine cylinders, rotational speed of turbo-blower, 
exhaust gas temperature before and behind the turbo-blower, 
supercharging air pressure and temperature, ship draught 
the integrated computer-based system for measuring and re- 
cording operational parameters of ship power plant devices 
- NORLIGHT Integrated Ship Monitoring and Control, I/O 
Station A 01 ER, NORIS Tachometer GMBH & Co. 

> for the measuring of ME fuel consumption — a calibrating 
tank of 11dm* capacity being an element of the standard 
equipment of ship power plant, and a stop-watch. 
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By making use of the standard measuring instrumentation 
its usefulness had to be assessed regarding its accessibility 
and easiness of read-out taking, and measurement accuracy of 
selected parameters. 

The accessibility and easiness of read-out taking is obvious. 
The standard measuring systems are fitted with easily readable 
indicators, usually in an analogue or digital form or in the form 
of messages displayed on a computer monitor. 

In the case of investigations of such complex technical ob- 
jects as ships, indication and measurement accuracy should be 
considered in a way different from that used in direct metrolo- 
gical measurements. In experimental investigations of complex 
objects the variability of measured quantity is assumed to be 
itself a basic source of errors [2]. The difference results from 
an influence of disturbing factors on the investigated object. 
Additionally, errors of the measuring method itself (including 
read-out errors) should by taken into consideration. In Fig.1 
is shown the example distribution of measurement results 
within a sample, location of the mean value of the sample, as 
well as the range of changes resulting from variability of the 
object and inaccuracy of the measuring method. In the case of 
long-lasting measurements (e.g. of fuel consumption or a large 
number of repetitions) one should also analyze errors resulting 
e.g. from a change of service conditions during measurement 
taking. The whole range of values of resultant factors is called 
the variability measure of results [2]. 

The errors resulting from the above mentioned sources can 
be numbered among random ones. As usually assumed, they 
are normally distributed that makes it possible to state that in 
accordance with the Gauss law the average from the sample is 
the most probable value of the resultant factor obtained from 
multiple measurements : 


(1) 


z, - a single measurement result 
r - number of repetitions of a measurement. 


Fig. 1. Diagram of position of variability measures 
and measurement error for an example sample of results : 
Zz, - a single measurement result, Òa - measurement method error, 
6. - variability measure of an object, 6, - variability measure of results, 
o ee zw 
Zj - average value from a sample. 


Possible systematic errors resulting from individual features 
of measuring instruments should be eliminated by introducing 
some corrections to read-out values on the basis of calibration of 
the measuring instruments. Also, should be rejected the results 
deemed to be loaded by rough errors, i.e. those significantly 
differing from the average from read-out values, resulting from 
mistakes or other hardly identifiable causes. 

The basic variability measure of measurement results is the- 
ir variance [2]. In the case in question, is of interest the variance 
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of the whole set of results coming from realization of the scope 
of experiment. For the same number of repetitions, r, for every 
point of the scope the following can be written [2, 3, 4] : 


o°(z) = = mira 7 EDALE- 


u=l j=l 


z'ye 2) 


where : 


m — number of measurement points 
r — number of repetitions in a measurement point 
z® — result of a single measurement 


Z® — average value in a measurement point. 


The root-mean-square error of results or standard 
deviation is equal to the square root of variance [2] : 


o(z) =./o°(z) = mye 


u=l j= 
The standard deviation of the average, 1.e. 
the mean standard error of the average 
is described by means of the following expression : 


Zz) 


o(z) 


o= = Se (u) zW)? 


ae Di = j=l 


(4) 


and the relative mean square error, i.e. the variability 
coefficient of results — by the expression : 


_ O(Z) 


~ (u) 
Z šr 


100% 


(5) 


where : 
me average from a sample of results. 


In the case of the direct measurements of the independent 
quantities x,, x,, ... x,, the root- mean-square error 
is equal to [5]: 


(2) =Jo(x) + 07K) +..40%(K,) © 
whereas the mean square error of the functions 
of direct measurements of the kind : z = X, X,,... X >) 
is equal to [5]: 
f f (2 
| ù | o°(x,) + | o°(x5)+... 
6(X,) 6(x>) 
(7) 
o(z) = sf 5 
: on o°(x,) 
5(x,) 


Discussion on usefulness of the standard measuring in- 
strumentation intended for energy measurements on ships was 
carried out on the example of measuring the ship speed and ME 
fuel consumption. Accuracy of measurements and measurement 
method was assumed an assessing criterion. 

The measurements were performed in accordance with 
a static, determined, multi-factor, rotational and uniform 
program of the technical experiment, described in [2], which 
contains 13 measurement points. The measurement accuracy 
was assessed in compliance with the above presented principle 
of interpretation and by using the formulae (2 + 16). 


ACCURACY ASSESSMENT 
OF SHIP SPEED MEASUREMENTS 


The accuracy assessment of ship speed measurements was 
performed for the results obtained from two measurement series 
(no. 2 and 4 according to the same notation as that used in me- 
asurement reports from shipboard tests), carried out in similar 
sea conditions, which are presented in the table below: 


No. of 
measurement 


Ship speed v[kn] 


Series no. 2 | Series no. 4 | On average 


After calculations by means of the formulae (2), (3), (4) 
and (5) the following results were obtained : 


variance of measured values of ship speed 07(v) = 0.0258 
standard deviation of o(v) = 0.16 knot 

standard deviation of the average value o; = 0.032 knot 
variability coefficient of measured values o = 1.8%. 


ACCURACY ASSESSMENT OF ME FUEL 
CONSUMPTION MEASUREMENTS 


de De be le 


The accuracy assessment of fuel consumption measure- 
ments was performed for the results obtained in average sea 
conditions, presented in the table below: 


No. of 
measurement 


Fuel consumption per hour B, [kg/h] 


Series no. 4 


Be =2 


Series no. 2 


nest 


On average 


BO = z® 


After calculations by means of the formulae (2), (3), (4) 
and (5) the following results were obtained : 


variance of measured values of fuel consumption 
o°(B,) = 4.33 

standard deviation of o(B,) = 2.1 kg/h 

standard deviation of the average value 0; = 0.41 kg/h 
variability coefficient of measured values o = 2.0%. 


ACCURACY ASSESSMENT 
OF THE MEASUREMENT METHOD 
FOR SHIP SPEED 


As the ship speed measurement is direct, its accuracy 
depends on accuracy of measuring instrument and external 
disturbances influencing the investigated object. The mean 
square error of measurement method is determined by using 
the relation (6). For a single measured quantity (ship speed v) 
it takes the form : 


VVV WV 


2 
o(v) = Jo (y) (8) 
hence it is equal to the square root 
from variance of the measured values : 
ia, = 
ov) = > D(i-¥) (9) 
i=l 


where : 


v, —result of a single read-out of ship speed 
v -—average speed in a measurement point 
r — number of repetition in a measurement point. 


To calculate the mean square error, were taken the measu- 
rement results at the ME rotational speed n = 830 rpm and 
propeller pitch H=80% for the sea conditions in centrum of 
measurement series 2 and 4, presented in the table below : 


8.9 
Number of repetitions r | 6 


After calculations by using the formulae (9) and (5) 
the following was obtained : 


On average V 


* the mean square error of the ship speed measurement 
method o(v) = 0.12 kn 
* variability coefficient of measurement results o = 1.3 %. 


ACCURACY ASSESSMENT 
OF THE MEASUREMENT METHOD 
FOR ME FUEL CONSUMPTION 


The ME hourly fuel consumption is the direct measurement 
function described by the formula : 


3600 V 
B,= - Pr 
where : 


[kg/h] (10) 


V = 11 dm’ — volume of calibrating tank 
P, = P,; — (t — 15)0.00064 [kg/dm’] — fuel density 
in temperature of measurements 
P,, = 0.875 kg/dm’ — fuel density at 15°C (standard density) 
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t [°C] — temperature of fuel 
t[s] — consumption time of fuel contained 
in calibrating tank volume. 


After inserting constant values the formula takes the form : 


35030.6 — 25.34t 
B, = z 


[kgh] (1) 

The mean square error of the fuel consumption measure- 
ment method is calculated on the basis of the relation (7). In 
the case in question it is described as follows : 


2 2 
o(B,,) = | S o°(t) + o°(t) 


Values of the partial derivatives of the function (11), with 
respect to fuel temperature and measurement time, after inser- 
ting the average measured values (see the table below), are as 
follows, respectively : 


= (12) 


OB, _ -25.34 _ 

e = m =- 0.07542 (13) 
OB,  25.34t-35030.6 
t= : =-0.3012 (14 
OT t 
The mean square error 
of fuel temperature measurements is equal to : 
1 < = 
o(t) = | — S(t, - 1)” (15) 
tele 


where : 


ti — result of a single read-out of fuel temperature 
t — average result of fuel temperature measurement 
r — number of repetitions in a measurement point. 


The mean square measurement error of consumption time 
of fuel contained in calibration tank volume is equal to : 


a= | A Sinn) 


where : 


(16) 


t, — result of a single read-out of consumption time of fuel 
contained in the calibrating tank volume 

T — average consumption time of fuel contained in calibrating 
tank volume 

r — number of repetitions in a measurement point. 


To determine the mean square error of fuel consumption the 
measurement results at the ME rotational speed n = 830 rpm 
and propeller pitch H = 80% for the sea conditions in centrum 
of 2 and 4 measurement series, presented in the table below, 

were applied : 


No. of 
measurement 
series 


On average 


Average time T [s] 


Number of repetitions r 
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After the calculations by using the formulae (12) and (5) 
the following was obtained : 


æ the mean square error of the fuel consumption measurement 
method o(B,) = 1.44 kg/h 
æ variability coefficient of results o = 1.4 %. 


ACCURACY ASSESSMENT 
OF MEASUREMENTS 
OF SHIP PROPELLER SHAFT TORQUE 


Torque measurements were carried out with the use of 
special instruments installed by the measuring team. For the 
measurements was used the PHILIPS PR 9914/01-NC 9408499 
strain-gauge torque meter with analogue read-out system. 

Accuracy assessment of torque measurements was per- 
formed for the results obtained from two measurement series 
(no. 2 and 4 - according to the same notation as that used in 
the measurement reports from shipboard tests) carried out in 
similar conditions, presented in the table below : 


No. of 
measurement 


Torque M[Nm] 


Series no. 2 | Series no. 4 | On average 


1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 


After calculations by using the formulae (2), (3), (4) and (5) 
the following was obtained : 


variance of measured torque values o7(M) = 18003.8 
standard deviation of o(M) = 134.2 Nm 

standard deviation of the average value ©; = 26.3 Nm 
variability coefficient of results o = 2.1%. 


USEFULNESS ASSESSMENT 
OF STANDARD INSTRUMENTS 
FOR ENERGY MEASUREMENTS 


HH HH 


O The obtained results of the assessment of measurement ac- 
curacy, performed with the use of ship standard instrumen- 
tation, are contained within the interval of 2% variability 
coefficient, and the assessment of accuracy of measurement 
methods — of 1.5%. It makes it possible to conclude that the 
standard measuring instruments are fully useful for carrying 
out energy measurements on ships. 


O The variability coefficient of measurement results obtained 
by means of the special instruments for torque measuring, 
amounts to 2.1%, hence it is not greater than that in the case 
of the standard measuring instruments installed on ships. 


Q 


the accuracy of the measurement methods, which is higher 
than that of measurement results, confirms that the reaso- 
ning presented in this paper is correct. 


Similar energy measurements were carried out in the labo- 


ratory of combustion engines of Ship Power Plant Department, 
Gdynia Maritime University. As expected, the obtained results 
were found even more accurate in view of the stable operatio- 
nal conditions and loading state of the laboratory engine. For 
instance the variability coefficient value of measurement results 


of 


hourly fuel consumption, equal to 0.65% was obtained. 
To conclude one should state that the standard measurement 


instrumentation which is installed as a rule onboard contempo- 
rary ships, can be successfully used for energy measurements. 


In 


particular, the measurements performed with the use of the 


standard instrumentation are of the following features : 


> 
> 


> 


> 


a sufficient accuracy 
no special instruments are required to be installed 

(a decrease of cost) 
they can be carried out by ship’s crew itself hence no ad- 
ditional measuring team is required (a decrease of cost) 
measurement results may be used in operational and dia- 
gnostic analyses of technical state of ship’s propulsion 
system. 


NOMENCLATURE 

B, [kg/h] - hourly fuel consumption 

H [%] - pitch of CP propeller (Controllable Pitch Propeller) 
KBN - The State Committee for Scientific Research 

M [Nm] - torque 

n [rpm] - rotational speed 

r - number of measurement repetitions 

t [°C] - temperature 

u - number of a measurement point 

v [kn] - ship speed 


V [dm?}] - volume 


- result of a single measurement 


Zj - average value from a sample of measurement results 

Zz - result of a single measurement point 

Zz - value in a measurement point 

Zo - average value in a measurement point 

ze - average value from a sample of measurement results 

Òn - measurement method error 

4, - variability measure of an object 

Ò - variability measure of results 

p [kg/dm?]- density 

(o - relative mean square error 

o - variance 

ofz) - mean square error 

o, - standard deviation of average value 

t [s] - time 
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Calculation of the mean long-term 
service speed of transport ship 


Continuation ot the first part of the paper published 
in the Polish Maritime Research No.4(50), October 2006 


Part ll 
Service speed of ship sailing 
on regular shipping route in real weather conditions 


Tadeusz Szelangiewicz 
Katarzyna Zelazny 
Szczecin University of Technology 


ABSTRACT 


Service speed obtainable by a ship in real weather conditions when sailing on a given shipping route, is 
one of the major parameters which have great impact on ship operation costs. The so far used, very ap- 
proximate method of service speed prediction based on “service margin”, is very little exact. In this paper 
a new method based on additional ship resistance dependent on mean statistical parameters of wave and 
wind occurring on a given shipping route, is presented. The mean long-term service speed is calculated on 
the basis of the calculated additional resistance and the screw propeller and propulsion engine parameters. 
Also, a new definition of service margin and a way of its calculation is presented apart from the results of 
the mean service speed calculation depending on ship 5 type and size and shipping route. 


Keywords : ship service speed, wind, waving, shipping route, service margin, long-term prediction. 


SHIP PROPULSION SYSTEM 
Screw propeller 


Ship propulsion system is represented first of all by a propel- 
ler cooperating with propulsion engine. To calculate instanta- 
neous ship service speed it is necessary to know hydrodynamic 
characteristics of the propeller (derived from model tests or 
calculated by using approximate methods) and the characteri- 
stics which describe engine load area. 

In the calculations were used the characteristics of B-Wa- 
geningen screw propellers [4,5], from which propeller’s thrust 
and torque were calculated. 

Thrust of the propeller placed behind the ship hull 
is expressed as follows : 


T= — 
It (27) 
where : 


Ro -~ total resistance of ship sailing in rough waters 
t — thrust deduction factor. 


Thrust of free propeller can be calculated from the formula : 
4.2 
T=K,p,D,1, (28) 
where : 


D_ — propeller diameter 
P : 
n, — propeller rotational speed 
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K, — thrust coefficient which — for typical B-Wageningen 
propellers of given values of the pitch ratio (P/D), 
blade area ratio (A,/A,), number of blades (Z) — is 
approximated by using the expression : 


K,=AptA J+ A,J7+ A? (29) 
where : 
Ap A, A, A, — coefficients of the polynomial describing 


thrust characteristics, depending on (P/D), (A,/A,), Z 
J — advance coefficient : 
J= ue w(V)I (30) 
pip 
w(V) — wake fraction depending on the ship speed V. 
The propeller — when working — generates the torque Q : 
= 542 
Q=KgpyD5n, (31) 
where : 
K,- torque coefficient which — like the thrust coefficient 
—can be expressed for a given propeller 
in the following form : 
— 2 3 
Ko = B, F BJ + BJ $ BJ (32) 
where : 


Bə B,, B,, B, — coefficients of the polynomial describing 
torque characteristics, depending on (P/D), (A,/A,), Z. 


Free-propeller efficiency (free from ship’s hull) is equal to : 


= Ar. J (33) 


The hydrodynamic characteristics of propeller achieved 
from model tests [4,5] are determined for a given Reynolds 
number and propeller surface state. In order to use the cha- 
racteristics for the behind-the-hull propeller were introduced 
relevant corrections associated with influence of Reynolds 
number and real state of propeller surface, which makes it 
possible — together with knowledge of ship hull surface state 
— to investigate influence of ship and propeller ageing on ship 
service speed. 

The expressions (27) and (28) for propeller thrust are valid 
for ship sailing in still and rough waters provided its oscillating 
and relative motions are so small that propeller emerging does 
not occur. During sailing in waves at large oscillating and rela- 
tive motions the propeller operates in highly aerated water or 
emerges. It generates thrust variations and drop of mean effec- 
tive thrust value relative to that in still water (even if the ship 
goes at constant speed and number of propeller rotations). 

The thrust decrease is caused a.o. by an influence of water 
particles being in oscillating motion, on wake, as well as by 
propeller emerging resulting from large relative motions of the 
ship in waves. The thrust decrease due to ship motions in waves 
has been discussed in various publications where approximate 
formulae for estimating the influence of relative motions on 
operational parameters of propeller are included [3]. 


The propeller thrust decrease due to ship relative motions in 

waves is defined by means of the following coefficient [3] : 

hp 

Kry [| 
Ky 


where : 


(34) 


Br= 


Ky (=) —thrust coefficient of emerging propeller (the 
quantities h, and R are shown in Fig. 11) 
K, ~thrust coefficient of fully immersed propeller. 


Changes of values of the coefficient p. 


h 
in function of (>) are shown in Fig.12, [6]. 


Fig. 11. The propeller draught h, as defined in the equation (34) . 


Thrust decrease due to propeller emerging 


Thrust decrease due to : 
- propeller emerging 

- stern wave system generating 

- variable lifting force on propeller blades 


0.5 4 


0.5 1.0 
Fig. 12. Thrust decrease during propeller emerging [6] . 
In simulative propeller thrust calculations, the influence 
of wave parameters on propeller operation effectiveness was 
taken into account (detailed analysis of the phenomenon is 
given in [9]) whereas the propeller emergence was mitigated 
by intentional reduction of ship speed. 


Propulsion engine 


The behind-the-hull propeller loads the ship engine 
with its torque (31). The relation of the propeller torque 
and propulsion engine output power is as follows : 


Pp 
= 7 35 
Q 2mm om 

where : 

P, ~ power delivered to the propeller 
n — rotational speed of engine (for slow-speed engines 
if no reduction gear is applied : n= n) 

Pa~ Nnw ie ie (36) 


N — propulsion engine output power 

Ng ~ relative rotative efficiency 

Nıw 7 shaft-line efficiency 

Np ~ reduction gear efficiency, if applicable. 


Engine output power for a designed ship is so selected as to 
obtain the propulsion system working point of the value equal 
to 0.85 N, at the design (contractual) speed in still water. In 
real weather conditions, when wind and waves affect the ship 
and the additional resistance AR appears, then the propulsion 
system working point changes its location within the propulsion 
engine’s load area. By controlling fuel charge (and this way also 
rotational speed of engine and propeller) the working point can 
be positioned in the maximum continuous rating area or the 
limited rating area for engine overloading (Fig.13). 

Therefore, in order to find an instantaneous speed at which 
a given ship would sail in considered weather conditions it is 
necessary to know the engine load area which is constrained by 
relevant characteristics. This area indicates where the propul- 
sion system working point is to be located. For instance the 
area of a Sulzer engine is shown in Fig.13. 

The particular areas are limited by the engine performance 
characteristic curves in the following form : 


N=k_-n™ 
m 
where : 
N — engine output power 
k „~ coefficient for a given characteristic curve 
n — engine rotational speed 


(37) 
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m — exponent depending on a type of engine of a given 
producer : 
- e.g. for SULZER RTA 52, RTA 62, RT 72, RTA 84 
slow speed engines it is equal to : 
m=0 -for constant rated or maximum output power 
m=1  -—for constant torque characteristic curve 
m = 2.45 — for overloading characteristic curve. 
Contract maximum continuous rating - CMCR 
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M, - nominal torque 
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Relative engine output power [% CMCR] 


®© - continuous rating area (marked green) 


- engine overloading area (red) 
© - engine sea trial operation area (deep blue) 


®© - engine optimum still-water operation area (blue) 
© - engine instantaneous operation area (yellow) 
Fig. 13. The load area of an example Sulzer propulsion engine [2] . 


Service speed of ships fitted with engines of other producers 
can be calculated if engine characteristics limiting 
its load area are known. 

The searching for of a working point of propulsion system 
and its parameters can be effected by taking into account vario- 
us criteria e.g. that for maintaining a given or maximum speed, 
or also by using the criterion for minimum fuel consumption 
or maximum efficiency of whole propulsion system. 


Intentional reduction of ship speed in view 
of hazardous wave-generated phenomena 


During ship sailing in rough water ship’s oscillating mo- 
tions and their derivatives i.e. velocities and accelerations are 
the direct effect of waves. The secondary phenomena which 
accompany the oscillating motions are : shipping of water on 
deck, propeller emerging, pounding of wave to ship bottom 
and sides (slamming), worsening of ship stability and ma- 
noeuvrability, additional dynamic loads on hull structure. The 
oscillating motions and accompanying phenomena depend on 
parameters of ship hull and waves as well as the ship speed V 
and course relative to waves (the angle B). The phenomena, 
especially if very intensive, may lead directly to averages and 
disasters at sea. It is possible to mitigate the phenomena, e.g. 
ship rolling - by changing ship’s course relative to waves (the 
angle p), reducing ship’s speed V, or simultaneous changing 
ship’s course and speed. 

In predicting the mean service speed was made the assump- 
tion that ship’s speed reduction will be performed if : 


Uz > Uzaop (38) 


where : 


Uz — mean statistical value of the wave-generated phe- 


nomenon Z considered hazardous to ship 
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Uzaop — permissible value of the wave—generated pheno- 
menon Z, at which the sailing ship is still safe. 


In the performed calculations the ship course relative to 
waves was kept unchanged (though its changing is always 
possible), since in predicting the mean long-term ship service 
speed it was assumed that the ship heads the course resulting 
from the selected shipping route. 

In assessing ship performance in waves and making deci- 
sion on possible reduction of its speed, the following pheno- 
mena were taken into consideration: rolling, pitching, vertical 
accelerations, horizontal transverse accelerations, shipping of 
water on deck, slamming, propeller emerging. The calculation 
methods of the phenomena and the assumed acceptance criteria 
are presented in [1, 5, 8]. 

The points at which values of the parameters of ship sea- 
-keeping qualities were calculated, are presented in Fig. 14, and 
in Fig.15, in the polar diagrams is shown an influence of ship 
speed and wave direction relative to ship on some phenomena 
induced by waves. The zone of unsafe intensity, i.e. where the 
sea-keeping criteria are exceeded , is marked red; yellow colour 
warns about a real hazard; green colour means the safe zone. 
It can be observed that reducing the speed one can mitigate 

intensity of the phenomena in question. 


= 


3 


37 
1) 


Fig. 14. The points at which ship sea-keeping parameters were calculated : 
1 — for slamming, 2 — for water shipping onto deck, and bow accelerations, 
3 — for propeller emerging, 4 — for accelerations at wheelhouse. 


Wave : 
H,=85m ; T,;=10.1s ; p=0° 


Water shipping on deck 
V=17w 
course angle = 40° 


Rolling 
V=15w 
course angle = 135° 


Vertical accelerations 
at wheelhouse 
V=17.5w 
course angle = 60° 


Horizontal transverse 
accelerations at wheelhouse 
V=10w 


course angle = 100° 


Fig. 15. Influence course angle and speed of K1 ship on its selected sea- 
keeping qualities. Green - safe operation of ship. Yellow - warning on a real 
hazard (lower value of a relevant criterion is exceeded). Red - ships safety 
is endangered (upper value of a relevant criterion is exceeded) . 


A METHOD FOR PREDICTING MEAN 
STATISTICAL SERVICE SPEED OF A SHIP 
SAILING ON A GIVEN SHIPPING ROUTE 


Instantaneous service speed of ship 


During ship motion in rough water, apart from still-water 
resistance also additional forces due to wind, waves and possi- 
ble surface sea current act on the ship. These actions generate, 
apart from an additional resistance, a transverse force and 
a moment tending to rotate the ship around its vertical axis. 
The transverse force results in ship drifting, and the moment 
— in ship course changing. Ship’s passive rudder is to be laid 
appropriately to keep ship course constant under action of the 
external rotating moment. Under the assumption that the ship 
course has to be maintained , from the solution of the following 
non-linear equations : 


YAV) + Y,(V) + RCV,B) + Y,(V,B,5,) = 0 
M,(V) + MCV) + M(V,B) + M,(V,B,5,) = 0 


together with the relevant equations describing additional 
resistance forces, (presented below), for given values of the 
ship speed V, course angle y, wind parameters (V,, y,), wave 
parameters (H,, T, u) and possible sea current parameters 
(Ve > Yo) the following quantities can be obtained : 


B — ship drift angle 

6, Z ship passive rudder angle 

AR — additional ship resistance due to wind, 
waves, current and passive rudder 

R, — total ship resistance to motion. 


Schematic diagram of the complete calculation algorithm 
of B, 6, and R, for given values of the ship motion parameters 
(V , y) , wind parameters (V, , y,), wave parameters (H,, T, p) 
and current parameters (V_., Yc), was presented in [9]. 

The instantaneous ship service speed in variable weather 
conditions is calculated from the solution of the set of non- 
linear equations in such a way as to obtain the propulsion sys- 
tem working point laying within the engine continuous rating 
zone when total ship resistance and propeller thrust , as well 
as propeller torque and driving engine torque become equal to 
each other, respectively . Making use of the equations (4), (5), 
(27 + 36) one obtains the set of two non-linear equations : 


(Ag+ A,J+A,J°+ A,J°) p,Din?- as 2 
(40) 
N 
(By+B,J+B,J?+ B,J’) -Lw _ 9 
2 TPD n 
where : 
J — advance coefficient (30) 
R, — afunction of total ship resistance, depending on the ship 


speed V, ship course angle y, wave parameters H,, T, y 
and wind parameters V, , y, 

driving engine output power determined by its charac- 
teristics valid within respective intervals of its rotation 
number n. 


N - 


To solve the set of equations (40) and determine instan- 
taneous ship service speed it is necessary to know the total 
ship resistance R, which depends not only on the statistical 
parameters of waves and wind occurring on a given shipping 
route but also on current ship speed, course and drift angle. 
Since all the quantities depend on random parameters of waves 
and wind the total ship resistance should be calculated for all 


statistical parameters of waves and wind occurring on a given 
shipping route. The calculation algorithm of instantaneous ship 
speed is presented in [9]. 


Mean statistical service speed of ship sailing 
on a given shipping route 


During the long-term sailing of ship on an assumed ship- 
ping route the additional resistance due to wave action will 
depend not only on a wave height (and period) but also on 
geographical directions: of wave, p , and ship course , y. Also 
the additional resistance due to wind will depend, apart from 
the wind speed V,, on the directions y, and y . It means that 
the additional resistance and also ship speed will depend on the 
values of the parameters of waves (H, , T, , p), wind (V, , y,) 
and ship motion (V, y), which may occur on a given shipping 
route within a long period of time. 

In the case in question the surface sea current is consi- 
dered to be a determinate phenomenon of the mean speed V, 
and direction angle y,. If a ship is assumed to enter a region 
of large-scale surface currents then the current action will be 
included into ship still-water resistance. 

Therefore the occurrence probability of the additional 
resistance AR of a given value, as well as the speed V which 
can be reached at occurrence of that additional resistance, 
depends on : 


> shipping route and probability of staying the ship in parti- 
cular sea areas 

> the statistical parameters of waves , (H,, T,, p), and wind, 
(V, > Ya) and probability of occurrence of the parameters 
in given sea areas 

> probability of occurrence of the ship motion parameters 
(V, y). 


The probability of being the ship in a given situation 
when sailing in waves on a given route, is as follows : 


Pit i (41) 
where : 

fi — probability of staying the ship in the sea area A 

fs = probability of staying the ship in the sea area A du- 
ring the season S 

f — probability of occurrence of the wave direction p in 
the sea area A during the season S 

Dee — probability of occurrence of the wave of the para- 
meters (H,, T,), propagating from the direction u 

sere o probability of the event that the ship moves with the 


speed V and heads the course w, respectively. 


In a similar way can be expressed the probability p, of being 
the ship in a given situation associated with wind state. In the 
calculations of additional resistance due to wind and waves it 
was assumed that wind speed and wave height are mutually cor- 
related, hence p = p,. As the event of being the ship in a given 
situation described by (41) will result in generating an addi- 
tional resistance and achieving a determined speed, hence : 


Pw ~ Pr Py (42) 
where : 
Pp ~ Partial occurrence probability (in given conditions) of 
additional resistance 
p, — partial occurrence probability (in given conditions) of 


instantaneous ship service speed. 


Values of the additional resistance due to wind, R 


that due to waves, Ryw, 


xa > and 
depend on random parameters of wind 
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and waves. Therefore the same values of R,, and R,.,, can oc- 
cur for different values of the parameters V,,y,,H,,T,, p, 
V, y. For each of the values of this way calculated additional 
resistance, a value of ship speed is calculated (the criteria con- 
cerning sea-keeping qualities are simultaneously examined to 
execute possible speed reduction in order not to violate them 
at given wave and wind conditions). 
The total probability P_,, of achieving the ship speed V at 
a given value of the additional resistance AR, is as follows : 


Ny Ng Ny Dyr Ay Ny 


A= Sl p=1 H,T=I V= w=l 


Pry = P, LV; (AR,)] (43) 


where : 
V(AR,) — instantaneous ship service speed in function 
of instantaneous additional ship resistance 
Hh ghey Tog My ly numbers of sea areas (crossed by 
a given ship), seasons, wave directions, wave parameters, 
ship speeds and courses, respectively. 


By calculating the distribution function of occurrence 
probability, f(V,), of the instantaneous ship speed V, the mean 
long-term ship service speed for a given shipping route, can 
be determined as follows : 


ny 
> Pry, Vi(AR = const) 
V= i=l 
-= 


(44) 


n, — number of intervals containing similar values 
of the instantaneous ship service speeds. 


On the basis of the presented calculation formula for the 
mean ship service speed, (44), were performed relevant cal- 
culations and analyses for example ships and shipping routes, 
whose results will be presented in the 3“ part of the paper. 


NOMENCLATURE 


i 
Ao 
Ap A, A, A,~ coefficients of the polynomial describing thrust 
characteristics 
, B,, B,, B,— coefficients of the polynomial describing torque 
characteristics 
— propeller diameter 
— probability of staying the ship in the sea area A 
— probability of occurrence of the wave of the parameters 
(H, , T,), propagating from the direction pu 
f; — probability of staying the ship in the sea area A during the 
season S 
fy» f, — probability of the event that the ship moves with the speed 
V and heads the course y, respectively 
E — probability of occurrence of a given wind direction 
H — significant wave height 


— propeller blade area ratio 


0: 
P 


B 
D 
f, 
fir 


yy 


J — advance coefficient 

K, — thrust coefficient 

Kœ  ~ thrust coefficient of emerging propeller 

Ko ~ torque coefficient 

ka — coefficient of a given characteristic curve of engine 


performance 
— propulsion engine output power 
— nominal output power of propulsion engine 
— rotational speed of engine 
— nominal rotational speed of engine 
— rotational speed of propeller 
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Ps — power delivered to propeller 

Py ~ combined probability of reaching a given value of 
instantaneous ship service speed at occurrence of a given 
value of instantaneous additional ship resistance 

P, — partial probability of staying the ship in a given situation 

Pr — partial occurrence probability of instantaneous additional 
resistance (in given conditions) 

P, — partial occurrence probability of instantaneous ship 
service speed (in given conditions) 

Q — propeller torque 

Ro — total ship resistance to motion 

Raa Ra M,, ~ mean wind-induced forces exerted to going ship 

© (R,,- additional ship resistance due to wind) 

Rs Rr ,M,, — passive rudder forces (R,, - additional ship 
resistance due to rudder) 

Ryy > Rw» M,w ~ mean wave-induced drift forces (R - additional 

' ship resistance due to waves) 

T — free-propeller thrust 

T, — mean characteristic wave period 

t — thrust deduction factor 

U, — mean statistical value of the wave-induced phenomenon Z 

V — ship speed 

V, — wind speed 

Ve — sea current speed 

V — mean statistical ship service speed 

w — wake fraction 

Z — number of propeller blades 

p — ship drift angle 

Ba — coefficient of propeller thrust decrease 

Ya — geographical direction of wind 

Yo — geographical direction of sea current 

AR ~— additional ship resistance due to weather conditions 

k — passive rudder angle 

No — free-propeller efficiency 

Np — reduction gear efficiency 

Ne — propeller rotative efficiency 

Nw  ~ Shaft-line efficiency 

A — rudder aspect ratio 

u — geographical direction of wave 

Py — water density 

y — geographical direction of ship course 
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